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Abstract 



In this paper, we study the thermal relaxation in the one-dimensional self-gravitating 
system, or the so-called sheet model. According to the standard argument, the thermal 
relaxation time of the system is around Nt c , where N is the number of sheets and t c 
is the crossing time. It has been claimed that the system does not reach the thermal 
equilibrium in this thermal relaxation timescale, and that it takes much longer time for 
■ the system to reach true thermal equilibrium. 

Q\ \ We demonstrate that this behavior is explained simply by the fact that the relaxation 

0^ ' time is long. The relaxation time of sheets with average binding energy is ~ 20Nt c , 

(— i \ and that of sheets with high energy can exceed 1000Nt c . Thus, one needs to take the 

O average over the relaxation timescale of high-energy sheets, if one wants to look at the 

^ \ thermal characteristic of these high energy sheets. 

ST 

1 Introduction 

The one-dimensional self-gravitating many-body system was originally discussed mainly as 
a simple toy model to understand the violent relaxation, Q| because the thermal relaxation 
timescale of its discrete realization, the sheet model, was believed to be long. Until 1980s, it 
had been generally accepted that the thermal relaxation time of the system of N equal-mass 
sheets is of the order of N 2 t C: where t c is the crossing time of the system. 

However, by means of numerical simulation Luwel et al.Q have demonstrated that the 
relaxation time is of the order of Nt c . Reidel and Miller 0, 0] reached a similar conclusion, 
though they reported the presence of systems which apparently did not relax for much 
longer timescale. 



In a series of papers, Tsuchiya et al. pi, 11, 12] have studied the thermal relaxation 



process of one-dimensional self-gravitating systems in detail, by means of the numerical 
integration over very long timescale (some of their experiments covered 5 x 10 8 t c ). They 
claimed that the thermal relaxation of the sheet model proceeds in a highly complex manner. 
In the "microscopic relaxation timescale" of Nt c , each sheet forgets its initial condition, and 
the system is well mixed. However, according to them, the system does not really reach the 
thermal equilibrium in this timescale, and the distribution function remains different from 
that of the isothermal state. They called this state a quasiequilibrium 



By pursuing the time integration for much longer timescale, Tsuchiya et al. [12] found 
that the system exhibits the transition from one quasiequilibrium to another, and they 
claimed that the thermal equilibrium is only realized by averaging over the timescale longer 
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than the timescale of these transitions. Thus, they argued that there exists the timescale 
for "macroscopic" relaxation, which is much longer than the usual thermal relaxation (what 
they called "microscopic relaxation"). 

In this paper, we try to examine the nature of this "macroscopic" relaxation of the 
one-dimensional sheet model. In section 2, we describe the numerical model. In section 
3, we present the result of the measurement of the relaxation time. It is shown that the 
relaxation time, defined as the timescale in which individual sheets change their energies, 
depends very strongly on the energy itself, and is very long for high energy sheets. This 
strong dependence of the relaxation timescale on the energy naturally explains the apparent 
"transient" phenomena observed by Tsuchiya et oZ.[12| Section 4 discusses the implication 
and relevance of our results. 



2 The Model 

2.1 Sheet model 

The Hamiltonian of the sheet model is given by 

N 

# = f E ^ 2 + 2 * Gm2 E N - *;l> (!) 

i=l i<j 

where X{ and Vi are the position and velocity of sheet i, m is the mass of the sheets, N is 
the number of the sheets and G is the gravitational constant. The crossing time is defined 

as 

1 [ae 

t c = \ — , (2) 

where M = mN is the total mass of the system. Following Tsuchiya et al. [|l~2j] and others, 
we use the system of units in which M = 4E = AttG = 1. In this system, t c = 1. 

A unique nature of the one-dimensional gravitational system is that there exists the 
thermal equilibrium, unlike its counterpart in three dimensions. Rybicki [|| obtained the 
distribution function 
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where e is the specific binding energy defined as 

e = y + *(x)-*(0). (4) 

Here, ^(x) is the specific potential energy. This distribution function satisfies the relation 

/ 3M \ 2 / 3x \ 

ex p \-^ £ ) = sech ■ ( 5 ) 
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We performed the time integration of the system with iV = 16, 32, 64, 128 and 256. For 
all systems, the initial condition is a water-bag with the aspect ratio x max /v max = 2.5. 

2.2 Numerical method 

The important character of the sheet model is that one can calculate the exact orbit of each 
sheet until two sheets cross each other. Thus, we can integrate the evolution of the system 
precisely (except for the round-off error). This may sound like a great advantage, compared 



to the systems in higher dimensions whose orbits can be calculated only numerically. Instead 
of numerically integrating the orbit of each sheet, we can calculate the exact orbit for any 
sheet, until it collides with the neighboring sheet. Thus, by arranging the pairs using heap, 
we can handle each collision in log A?" calculation cost. 

Note, however, that typically each sheet collides with all other sheets in one crossing 
time. Thus, the calculation cost is 0(N 2 log N) per crossing time. Our simulation with 
N = 64 for 2 x 10 7 t c took 8 hours on a VT-Alpha workstation with DEC Alpha 21164A 
CPU running at 533 MHz. For this run, the total energy of the system was conserved better 
than 3 x 1(T 12 . 



3 Results 

3.1 Approach to the thermal equilibrium 

Figure [l] shows the time-averaged energy distribution function N(e), for different time 
periods and number of sheets. In all figures, the thin solid curve is the energy distribution 
of the isothermal distribution function of equation ([|). What we see is quite clear. As we 
make the time interval longer, the time- averaged distribution function approaches to the 
isothermal distribution. Thus, the numerical result suggests the system is ergodic. However, 
it also shows that the time needed to populate the high-energy region is very long. The 
sampling time interval is 128 time units for A" = 16, and 512 time units for A^ = 64 and 128. 
Thus, in the case of T = 2 18 and N = 16 (dash-dotted curve in figure la), total number of 
sample points is 2 15 = 32768. 

If we can assume that the sample points are uncorrelated, the possibility that no sample 
exceeds energy level £q is given simply by 

P(e ,N) = [l-P(s<s )] n , (6) 

where 

P(e <e )= / N(e)de, (7) 
J o 

and n is the number of sample points. Figure ^ shows 1 — P(e) as a function of e. For 
e = 1.25, P(s) = 0.996, and therefore the probability that none of 32768 samples does not 
exceed e = 1.25 is practically zero (< e~ 100 ). In other words, the numerical result seems 
to suggest that the system is not in the thermally relaxed state even after 2 x 10 5 crossing 
times. 

Of course, this result is not surprising if the relaxation time is long. Samples taken with 
the time interval shorter than the relaxation time have a strong correlation, and therefore 
the effective number of freedom can be smaller than n. Roughly speaking, if the relaxation 
time is longer than 10 , our numerical result is consistent with the assumption that the 
system is in the thermal equilibrium. In the next subsection, we investigate the relaxation 
time itself. 



3.2 Relaxation timescale 

We measured the following quantities: 



< £i(t ) - £i(t + At) > 

Dl = At ' (8) 

< h(t ) - e l {t + At)] 2 > 

D 2 = - . (9) 
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Figure 1: The time-averaged distribution function in the energy space N(E); (a) N = 16, 
(b) N = 64, (c) AT = 128. 




Figure 2: The compliment of the cumulative distribution function 1 — -P(e) for the thermal 
equilibrium. 

These quantities correspond to the coefficients of the first and second-order terms in the 
Fokker-Planck equation for the distribution function, and have been used as the measure 
of the relaxation in many studies (see, e.g., Hernquist and Barnes, ]]]] Hernquist et al. Q), 
for three-dimensional systems. However, to our knowledge this measure has not been used 
for the study of the sheet model. 

In order to see the dependence of these diffusion coefficients on the energy, we calculated 
them for intervals of Ae = 0.15. Figure ^ shows the results, for N = 16,64 and 256. The 
time interval At was taken equal to Nt c . We used smaller values for At and confirmed 
that the choice of At has negligible effect if At is larger than 10t c and smaller than 4iVt c . 
Time average is taken over the whole simulation period. We can see that both the first- 
and second-order terms show very strong dependence on the energy of the sheets, and of 
the order of l/100iV for e ~ 1.5. Figure || suggests that the relaxation timescale grows 
exponentially as energy grows. This behavior is independent of the value of N. 

We can define the relaxation timescale as 

t r = e 2 /D 2 , (10) 

that is, the timescale in which energy changes significantly. Figure |I| shows this relaxation 
timescale for different values of N and e. The relaxation time shows very strong dependence 
on the energy and the relaxation of high-energy sheets is much slower than that of sheets 
in lower energies. This is partly because of the dependence of t r on e itself. However, as w 
can see in figure ||, the dependence of the diffusion coefficient is the main reason. 

This result resolves the apparent contradiction between the fact that the relaxation 
timescale is of the order of Nt c H and that the system reaches the true thermal equi- 



librium only in much longer timescale. [12] It is true that the relaxation timescale is O(N), 
but the coefficient before N is quite large, in particular for sheets with high energies. 

An important question is why the relaxation timescale depends so strongly on the energy. 
This is provably due to the fact that high-energy sheets have the orbital period significantly 
longer than the crossing time. Typical sheets have the period comparable to the crossing 
time, and therefore they are in strong resonance with each other. However, a high-energy 
sheet has the period longer than the crossing time, and thus it is out of resonance with the 




Figure 3: The diffusion coefficients (a) D\ and (b) D2 piotted against the energy e for three 
values of N. Long-dashed, solid, and short-dashed curves are the results for N = 16, 64 
and 256, respectively. 




Figure 4: The relaxation time in unit of Nt c plotted against the energy e for three values 
of N. Curves have the same meanings as in figure [3| 



rest of the system. Therefore, the coupling between high energy sheets and the rest of the 
system is much weaker than the coupling between sheets with average energy. This explains 
why the relaxation of high energy sheets is slow. 



4 Summary and Discussion 

In this paper, we studied the thermal relaxation process in one-dimensional self-gravitating 



systems. We confirmed the result obtained by Tsuchiya et al.[V2\ that the thermal relaxation 
takes place in the timescale much longer than Nt c . However, we found that this is simply 
because the thermal relaxation timescale is much longer than Nt c . Even for typical sheets, 
the relaxation timescale is around 10Nt c . In order to obtain good statistics, we need to 
take average over many relaxation times. Moreover, the relaxation time for sheets in the 
high-energy end of the distribution function is even longer, since the relaxation timescale 
grows exponentially as the energy grows. Thus, it is not surprising that we have to wait for 
more than 10 4 iVt c to obtain good statistics. 

Does this finding have any theoretical/practical relevance? Theoretically, there is nothing 
new in our result. What we found is simply that numerical simulation should cover the 
period much longer than the relaxation timescale to obtain statistical properties of the 
system, and that the relaxation timescale of a sheet depends on its energy. Both are 
obvious, but some of the previous studies neglected one or both of the above, and claimed 
to have found a complex behavior, which, in our view, is just a random walk. 

Our finding of the long relaxation time by itself has rather little astrophysical significance, 
since in the large N limit, the relaxation time is infinite anyway. However, since any 
numerical simulation suffers some form of numerical relaxation, it is rather important to 
understand how the relaxation effect changes the system. To illustrate this, we examine the 



claims by Tsuchiya et al. [12] in some detail here. 

They argued that the evolution of the mass sheet model proceeds in the following four 
steps: (1) viliarization, (2) dynamical equilibrium, (3) quasiequilibrium, and (4) thermal 
equilibrium. According to them, the viliarization timescale is order of t c , and the energy 
of each sheet is "conserved" in the dynamical equilibrium phase, which continues up to 
t ~ Nt c . Then, "microscopic relaxation" takes place in the timescale of t ~ Nt c , where the 
energy of each sheet is relaxed, but the whole system needs timescale much longer to reach 
the true equilibrium, because of some complex structure in the phase space. 



Our numerical results are in good agreement with those of Tsuchiya et al, [12] but our 
interpretation is much simpler: First system virializes, and then relaxation proceeds in 
the timescale of thermal relaxation, which depends on the energy of the individual sheets. 
Thus, the central region with short relaxation time relaxes to the distribution close to the 
thermal relaxation in less than 100Nt c , but the distribution in the high-energy tail takes 
much longer to settle. In addition, the small number statistics in the high-energy region 
makes it necessary to average over many relaxation times to obtain good statistics. In other 
words, there are no distinction between the "microscopic" and "macroscopic" relaxation, 
and the evolution of the system is perfectly understood in terms of the standard thermal 
relaxation. 



Acknowledgments 

I would like to thank Yoko Funato, Toshiyuki Fukushige and Daiichiro Sugimoto for stim- 
ulating discussions, and Shunsuke Hozumi for comments on the manuscript. This work 
is supported in part by the Research for the Future Program of Japan Society for the 
Promotion of Science (JSPS-RFTP97P01102). 



References 

[1] L. Hernquist and J. E. Barnes, Astrophys. J. 349 (1990), 562. 

[2] L. Hernquist, P. Hut and J. Makino, Astrophys. J. Lett. 411 (1993), 248. 

[3] M. Luwel, G. Severne, and P. J. Rousseeuw, Astrophys. Space Sci. 100 (1984), 261. 

[4] D. Lynden-Bell, Mon. Not. Roy. Astron. Soc 136 (1967), 101. 

[5] G. Severne, M. Luwel, and P. J. Rousseeuw, Astron. Astrophys. 138 (1984), 365. 

[6] C. J. Reidl, Jr. and B. N. Miller, Astrophys. J. 318 (1987), 248. 

[7] C. J. Reidl, Jr. and B. N. Miller, Astrophys. J. 332, (1988) 619. 

[8] G. B. Rybicki, Astrophys. Space Sci. 14, (1971) 56. 

[9] L. Spitzer and M. H. Hart, Astrophys. J. 164, (1971) 399. 
[10] T. Tsuchiya, T. Konishi, and N. Gouda, Phys. Rev. E50 (1994), 2607. 
[11] T. Tsuchiya, N. Gouda, and T. Konishi, Phys. Rev. E53 (1996), 2210. 
[12] T. Tsuchiya, N. Gouda, and T. Konishi, Ap. SS. 257 (1997), 319. 



